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Abstract 

The critical behaviour of anisotropic Heisenberg models with two 
kinds of antiferromagnetically exchange-coupled centers are studied 
numerically by using finite-size calculations and conformal invariance. 
These models exhibit the interesting property of ferrimagnetism in- 
stead of antiferromagnetism. Most of our results are centered in the 
mixed Heisenberg chain where we have at even (odd) sites a spin-S (5") 
SU(2) operator interacting with a XXZ like interaction (anisotropy 
A). Our results indicate universal properties for all these chains. The 
whole phase, 1 > A > — 1, where the models change from ferromag- 
netic (A = 1) to ferrimagnetic (A = —1) behaviour is critical. Along 
this phase the critical fluctuations are ruled by a c = 1 conformal 
field theory of Gaussian type. The conformal dimensions and critical 
exponents, along this phase, are calculated by studying these models 
with several boundary conditions. 

1 Introduction 

The critical properties of one-dimensional regular Heisenberg spin chains 
with one kind of antiferromagnetically exchange-coupled spins have been 
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extensively studied in the literature. The prototype of these models is the 
anisotropic S = \ Heisenberg model or XXZ chain]!]]. This model is exactly 
integrable with a critical line of continuously varying critical exponents as 
we change the anisotropy (A), bringing the model from the ferromagnetic 
(A = 1) to the antiferromagnetic (A = —1) isotropic points. With the 
advance of the conformal invariance ideas the whole operator content of 
this model was obtained^, £|. The critical fluctuations are governed by a 
Gaussian type conformal field theory with conformal anomaly c = 1 and, 
moreover, the underlying currents satisfying a U(l) Kac- Moody algebra||. 

The extension of the XXZ chain to higher spins S > | attracted consid- 
erable attention after Haldane|| conjectured that, for the isotropic antifer- 
romagnetic point (A = —1), the model is critical or not depending if S is 
half-odd-integer or integer, respectively. Consistent with this conjecture, nu- 
merical calculations ]7[ |8|, [1 indicate that in the case of half-odd-integer spins 
the models are critical in the whole range of anisotropics (1 > A > — 1) 
from the ferromagnetic to the antiferromagnetic point. In the case where 
S is integer, a critical line starting at the ferromagnetic point ends at A c 
before the antiferromagnetic phase (A c > —1) entering the massive Hal- 
dane phase [0, §, |[|. For all spins the massless phases are ruled by a c = 1 
Gaussian-like conformal field theory [Q. 

In this paper we extend these studies by studying a quantum chain in 
which two types of antiferromagnetically exchange-coupled spins S and S' 
are located at alternate sites. When S ^ S', according to a theorem due 
to Lieb and Mattis[10], the isotropic model (A = —1) exhibits ferrimagnetic 
order, with a (S — S') -|-degenerate ground state, where L is the chain length. 
Consequently as we vary the anisotropy the model goes from the ferromag- 
netic point (A = 1) to the ferrimagnetic point (A = —1). We studied these 
models for (S,S f ) = (|, 1) and (S, S') = (|,|), using finite-size scaling and 
conformal invariance 0. To supplement our studies we considered two other 
anisotropic models that also exhibit ferrimagnetic behaviour at the isotropic 
point. Our studies show that all these models between the two (ferromagnetic 
and ferrimagnetic) isotropic points 1 > A > — 1 have a universal Gaussian 
critical behaviour with central charge c = 1. In this massless phase the 
critical exponents exhibit a model-dependent variation with the anisotropy. 

The paper is organized as follows. In section 2 we define the 
(S, S')— Heisenberg chain and review some results obtained from conformal 
invariance in the model with S = S'. In section 3 and 4 we present our nu- 
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merical results for the (S, S')— Heisenberg model and related models. Finally 
in section 5 we present our general conclusions. 



2 The model and conformal invariance rela- 
tions 

The mixed Heisenberg quantum chains are defined by attaching an SU(2) 
spin-S 1 at the odd sites (&i = (of , of , of ); i = 1,3,5,...) and a spin S' at 
the even sites (Si = (Sf, Sf, S'f); i = 2, 4, 6, . . .). The Hamiltonian on an L 
(even) site chain, with periodic ends, is defined by 

L/2 

H = - £ ( <-i s * + <-i<^ + &<-iS Z 2i ) (1) 
1=1 

where A is the anisotropy constant. This Hamiltonian, like the standard 
spin-5* XXZ chain (S = S'), has a U(l) symmetry due to its commutation 
with the z component of the total spin 

L/2 

S* = £ ("L-i + S*) (2) 

For A > 1 the model is massive and ferromagnetic with a double degenerate 
ground state corresponding to the two fully ordered states with S z = ±|f( S + 
S'). At A = 1 the lowest energy in all U(l) sectors 
S z = -f ( S + S' ), -| ( S + S' ) + 1, . . . , § ( S + S' ) - 1, f ( 5 + S" ) ) are de- 
generate, rendering a ferromagnet ground state with total spin ^(S + S') 
and a massless spectra with a quadratic dispersion relation. For 1 > A > — 1 
the ground state is single or double degenerate, depending if \S — S'\\ 
is integer or half-odd-integer and belongs to the sectors with S z = or 
S z = ±2, respectively. At A = — 1 the lowest energies in the sectors where 
S z = -\ \3 - S'\, -\ \S -S'\ + 1,...,%\S- S'\ -1,%\S- S'\ become de- 



generate and we have ferrimagnetic order |TIJ. For A < — 1 the ground state 
is double degenerate, occuring in the sectors with S z = ±^\S — S'\, and we 
expect a massive behaviour as in the standard S = S' XXZ chain. In order to 
illustrate the spectral dependence on the anisotropy A, in Fig. 1 we draw in 
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schematic form the location of the lowest eigenenergies of (|1|) in the various 
S z sectors. 

Our analysis indicates that in the whole region 1 > A > — 1 the model is 
critical, like the S = \ XXZ model (S = S 1 = ~). We assume that the Hamil- 
tonian ([!]), like most statistical mechanics quantum chains are conformally 
invariant in its critical regime. Under this assumption the machinery arising 
from conformal invariance tells us that for each primary operator 0, [Tl[ O a 
with dimension x a and spin s a in the Virasoro operator algebra of the infinite 
system, there exists an infinite tower of states, in the quantum Hamiltonian, 
for a periodic chain of L sites, whose energy and momentum as L — > oo are 
given by 

El y = E (L) + T v(x a + j+f) + o (iT 1 ) (3) 

and 

2-7T 

P-f = T (Sa+j-j') (4) 

where = 0, 1, — Here E (L) is the ground-state energy and v is the 
velocity of sound, which can be determined by the energy-momentum dis- 
persion relation or from the difference among consecutive energy levels in a 
same conformal tower. The finite-size corrections of the ground-state energy 
also give a way to calculate the conformal anomaly. For periodic chains, the 
ground-state energy behaves asymptotically asfl2| 



EoW -~ nC °+o(L-*) (5) 



L °° 6L 2 

where is the ground-state energy per site in the bulk limit. 

In the case where S = S' a critical phase appears |7|, |], |[ in (|l[) for 
anisotropies 1 > A > A c (5), where due to the Haldane conjecture^], A c = 
— 1 or A c > —1 depending if S is half-odd- integer or not. This massless 
phase is described by a U(l) conformal field theory with central charge c = 
10. The anomalous dimensions x a appearing through (0) in the U(l) sector 
S z = n of the Hamiltonian (ffl) with periodic ends are given by 

,2 



n 2 x 



m 



+ — , n,m = 0,±l,±2,... (6) 



4:Xp 



where x p depend on S and A. For S = S' = \ we have the exact dependence 
U Xp — (tt — cos _1 (— A)) /2tt. Although the model is not integrable for 
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S = S' > ~, numerical calculations indicate the conjecture || 
x p = (n — cos _1 (— A)) /4:7iS for — 1 < A < 0. Beyond the dimensions @ 
other integer dimensions also appear in the sector with S z = n = 0. This 
fact indicate that the underlying conformal field theory satisfies a larger 
algebra than the Virasoso conformal algebra, namely, a U(l) Kac-Moody 
algebra^, j8j. The dimensions @ correspond to operators O n m and the num- 
ber of its descendants will be given by the product of two U(l) Kac-Moody 
characters. The dimensions (0) indicate that the operators O n , m correspond 



to the Gaussian model operators |13| composed of a spin-wave excitation of 
index n and a "vortex" excitation of vorticity m. 

Other interesting properties of these c = 1 critical phases appear when 
we consider these chains with more general boundary conditions, compati- 
ble with its U(l) symmetry, i.e., by preserving the total spin S z as a good 
quantum number. Two of these conditions are, the x — y twisted boundary 
conditions 

Sl +1 ± iS y L+1 = (Sf ± iSf) (7) 

£>L+1 = -^lJ ( 8 ) 

where $ is an arbitrary angle, and free boundary conditions 

Sl +1 = Sl +1 = S Z L+1 = 0. (9) 

The net effect of the boundary angle $ in the dimensions (^|) is to shift 
the spin- wave index by an amount [5| $/27r, 

2 (m + $/27r) 2 . . 

Xn,m+$/2n = U X p H . (10) 

^±Xp 

In a semi-infinite lattice the correlation functions involving lattice points near 
the surface have a power-law decay distinct from the case where the points 
are away from the surface (bulk behaviour). These correlations are ruled by 
the surface exponents x s . These exponents can be obtained from the finite- 
size corrections of the mass-gap amplitudes of finite chains with free ends. 
Instead of (Q) and (f|), to each surface exponent of the semi-infinite system, 
at the critical point, there exists a set of states with energies given by[[Ll[| 

E<P = Ejp{L) + ^{x.+r) + o(L- 1 ) (11) 
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where Eq (L) is the ground-state energy of the L-site chain and r = 0, 1, 2, 
Instead of (||) we have|l2[ 



where foo is the bulk limit of the surface energy. The study of (jl|) with S = S' 
and free ends0, [14] shows that in the critical region, for each sector S z = n 
there appears only one conformal tower associated to the dimensions 

x s {n) = 2n 2 x p , n — 0,1,2,... (13) 

with the multiplicity of its descendants given by the character of a single 
U(l) Kac-Moody algebra g ITS . 



3 Results for the mixed Heisenberg chains 

We calculate numerically the eigenspectra of the Hamiltonian (|1|) by using 
the Lanczos method in the cases where (S, S') = (|, 1) and (S, S') = (|, |) 
up to lattice sizes L = 20 and L = 16, respectively. Our results, for several 
boundary conditions, indicate that in the whole region 1 > A > — 1 the model 
is gapless and conformally invariant. The ground state in this region will have 
the lowest possible value of \S Z \. It is non degenerate with S z = if — S'\ 
is integer and is doubled degenerate with S z = ±|, otherwise. Consequently 
in order to obtain a uniform convergence of our finite-size results we consider 
in the case (S, S') = (|, 1) only lattice sizes which are multiples of 4. 

Let us consider initially the case of periodic chains. The model is invariant 
under translation by a unit cell with two spins, with momentum p = {I = 
0, 1, . . . , |f — 1). In order to calculate the conformal anomaly and exponents 
from (0-0) we should estimate the sound velocity. As in the S — S' — | 
case the lowest eigenenergy with nonzero momentum ( modulo n ) belonging 
to the ground-state sector is associated to a primary spin-1 operator with 
dimension equal to unity for all values of A. Using Eq. (|3|) we obtain an 
estimate for the sound velocity 

(E^/L - E (L)) L 
v{L) = ± '— (14) 

Z7T 
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where E (L) is the ground-state energy and E^/l is the lowest eigenenergy 



of a state with momentum Air/L (modir). Using (14) the conformal anomaly 
c is obtained by extrapolating the numerical sequence obtained from (^). In 
table 1 we show, for some values of A, our estimates for c in the two cases 
{S,S') = (|,1) and (S,S f ) = (|, |). All the extrapolated results reported 
in this table, as in the subsequent ones, are calculated by using the alter- 



nating e-algorithm||16||, which is a variant of the van den Broeck-Schwartz 
method |17|. The errors are roughly estimated from the region of stability of 
these approximants. It was not possible to obtain reliable results near the 
isotropic points A = 1 and A = — 1. This happens because like near the 
ferromagnetic models with S = 5"[|J, the sound velocity decreases towards 
zero as we tend towards the isotropic point (the energy-momentum disper- 
sion relation changes from linear to quadratic). Our results indicate that we 
have a conformal anomaly c = 1 in both cases and we believe that this is the 
general case for arbitrary S ^ S', since the spectrum (see figure 1) shows the 
same essential features independently of S and S' being integer or half-odd- 
integer. Moreover the vanishing of the sound velocity as we tend toward the 
ferromagnetic (A = 1) and ferrimagnetic (A = —1) points indicate that the 
critical fluctuations around the ferrimagnetic point are similar to those near 
the ferromagnetic point. As we see in table 1 the results for (S, S') = (|, |) 
are slightly better than those of (S, S') = (|, 1). This is due to the fact that 
for (S, S 1 ) = (|, |) the number of terms in the finite-size sequences are larger 
since L can be an arbitrary even number. 

The conformal dimensions are calculated by using ([3]) and (|T4]). For 
example in the (S, S') chain the lowest energy E n in the sector with total 
spin S z = n is associated to the dimension x^ s , which is calculated from the 
asymptotic (L — > oo) value of the sequence 

T s,S' (n _ E n (L)-E Q (L) 

where v(L) is given by (p!4[). In tables 2 and 3 we show our results for n = 1 
and 2 and some values of the anisotropy. We see from these tables that for 
all values of the anisotropy, the extended relation 

x s n s ' = n 2 x S { S ' (16) 

holds. These dimensions are similar to the Gaussian dimensions x n fl appear- 
ing in (|J), on identifying x p = xf' s . Different from the critical regime in the 
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homogeneous spin case S = S', when S ^ S' the dimensions xf' 5 increase as 
we depart from the ferromagnetic point, but around A < — 0.5 it starts to 
decrease again and we have small values of the ferrimagnetic point, 

as we normally see near the ferromagnetic point. The small values of the 
exponents x t ' near the ferromagnetic and ferrimagnetic points is the sig- 
nature of the long-range ordered ground state at the isotropic points. The 
fourth column of these tables also shows that near the ferromagnetic point 
Xi' S = x}^ 2 ' 1 ^ 2 / (S + S'), which give us 

s>s , _ n - cos(-A) 
Xl - 2n(S + S') ' A (17) 

This result when compared with the conjectured || results for S = S', 
indicate that near the ferromagnetic point we have essentially a Heisenberg 
model with effective spin (5* + S')/2. On the other hand the degenerascy of 
the ground state at the ferrimagnetic point (A = —1) induce us to expect 
near this point an effective Heisenberg ferromagnetic chain with effective spin 
(S — S')/2. However, we are not able to make a conjecture like in Eq. (|17]). 

1/2 1 

We also made an independent calculation of the exponents x± ' by us- 



ing standard finite-size scaling|18|. This exponent is related to the ratio of 
the "electric" susceptibility 7 P |19| and correlation-lenght exponent z/, by the 
relation 

^ = 2(l-x\ /2 ' 1 ). (18) 

The "electric" susceptibility is the response of the system to an staggered 
transversal field. This susceptibility Xl * s calculated by adding an "electric 
field" interaction \9A^i + QbS^ in (|) 



(19) 



where -Eo(0 is the ground-state energy in the presence of the "electric" field. 
The Lande factors qa and gs produces the staggering effect of the transverse 
"electric" field. In the isotropic case S = S' we must choose Qa/qb ^ 1, 
otherwise it will produce the effect of a uniform transverse field, which is 
not related with the exponent j p . In the case S ^ S' our results shows that 
a similar effect also occur and in order to calculate 7 P we should consider 
Qa/qb 7^ r c . For S = 1/2 and S' = 1, r c changes from r c = 2 to r c = 
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2.66 as the anisotropy changes from the ferromagnetic point (A = 1) to 
the ferrimagnetic one (A = —1). These points are probably related to the 
compensation mechanism which usually happens in ferrimagnetic ordered 
models when temperature effects are taken into account [2(J. In the bulk 



limit Xl ~ L lp l v . The extrapolation of the finite-size sequence obtained by 
choosing = and gs = 1 gives from (|18|) the results in the last column of 
table 2. In this case, since we calculate lattices up to L = 16, it is difficult 
to obtain an error estimate through the alternating e-algorithm [|l6)| . The 



results presented are in reasonable agreement with those derived by using the 
conformal invariance relations. Beyond the dimensions presented in tables 2 
and 3 our results also indicate other dimensions which would correspond to 
x n , m in (|6|) with m/0. Instead of presenting these dimensions we show in 
table 4 the lowest dimensions x$ s obtained by calculating the (S, S') model 
with the twisted boundary conditions given by (|7||f). These dimensions are 
obtained from the bulk limit extrapolations of the sequence 

s,s> (n _ E*(L)-E (L) 

where E$(L) is the ground-state energy of the Hamiltonian (P with L sites 
and boundary conditions (HB). The results given in table 4 indicate that 

x s,s> _ (^) 2 , 21) 
~ 4xf' s ' 

where xf' 5 is given in tables 2 and 3. These dimensions correspond to the 



dimensions Xq^/2w in ([TOf) . The results presented in tables 2-4 clearly indicate 
that in the whole disordered regime — 1 < A < 1 the conformal dimensions 
are those of a Gaussian model. The dimensions are given by © where 
Xp = X\ is a continuous function of A with some of its values given in 
tables 2 and 3. We have also studied the (S, S') = (|, 1) model with lattice 
sizes L = 41 + 2 (I = 1 — 4). In this case the ground state is degenerate like 
in the standard 5" = S' — | XXZ chain with an odd number of sites [0] and 
we obtain the dimensions x n+1 / 2 .m ( n , rn — 0, ±1, ±2, . . .). 

For completeness, we also calculated the surface exponents of these chains. 
These exponents are calculated from the eigenspectra of (|l|) with free bound- 
ary conditons. From (|TT|) the surface exponents xf ,s (n) associated with the 
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state with lowest energy E„ q in the sector S z = n of the (S, S') chain can 
be estimated from the large- L behaviour of the sequence 

F (n,L)--^ (22) 

-^0,1 — -^0,0 

where E^ is the m excited state in the sector S z = n. The estimator (p2]) 
was obtained by assuming that like in the S = S' case, the first mass gap 
amplitude in the ground-state sector is associated with a dimension equal to 
unity. Our results are shown in table 5 where we clearly see the same type 
of extended relation as in (|16D, namely 

xf s '{n) = n 2 xf s '(l). (23) 

Comparing these results with those of tables 2 and 3 we obtain the relations 
(|13"D expected in a Gaussian model 

xf s '(n) = 2n 2 x s { s ' , n = 0, 1, 2, . . . (24) 

where xf' s is the dimension which appeared in the periodic case. 



4 Results for other related models of ferri- 
magnetism 

Inspired by the results of last section we will try to see if the general critical 
features of the (S, S^-Heisenberg models can also be observed in other models 
exhibiting ferrimagnetism. In this direction we will study two other models 
defined on a bipartite lattice as in figures 2a and 2b. At each lattice point 
we attach a spin-| operator which interacts along the lines of Fig. 2 with 
interactions of XXZ type 



H = - E ( of at + afal + AofoJ ) . (25) 

<ik> 

The first model defined in the bipartite lattice of figure 2a we denote ABC 
and the second one in figure 2b we denote AB2. The lattice size L is consid- 
ered as two times the number of lattice sites in the sublattice A (in figures 
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2a and 2b L = 6). Both models, at the isotropic ferromagnetic point A = 1 
are fully ordered. At A = — 1 they show a ferrimagnetic behaviour since the 
total number of spin variables in each sublattice is not equal. As before the 
Hamiltonian has a U(l) symmetry and its Hilbert space is separated in the 
<7 2 -basis into block disjoint sectors labelled by the ^-component of the total 
spin S z = J2i<?i- The ground state location in these sectors as well as its 
degeneracies on a finite lattice for — 1 < A < 1 are exactly like those shown 
in figure 1, on taking S = \ and S' — 1. 

Our numerical results for periodic boundary conditions indicate that both 
models are disordered and massless in the whole regime 1 > A > — 1. The 
thermal effects of the model AB2 at A = — 1 are considered in [^T|. We now 
report our results separately for both models. 



4.1 Model ABC 

Using equation ( |T4|) we observe as in section 3 that the sound velocity ap- 
proaches zero as we tend towards the isotropic points A = ±1. The conformal 
dimensions x^ BC (n = 1, 2, . . .) associated to the lowest eigenenergy in the 
sector S z = n are obtained by extrapolating the sequence (|T5| ) for L up to 
16. Our results for some dimensions and anisotropies are shown in table 6. 
As we see in this table, like in (|T6| ) the relation x^ BC = r?x^ BC also holds 
for the whole range of anisotropies 1 > A > — 1, with x^ BC depending con- 
tinuously on A. Since at A = 1 we have the same ground-state degeneracy 
as in the (S, S') = (~, 1) Heisenberg model we would expect an asymptotic 
behaviour, as A — > 1, like ([[7]) with S = | and 5" = 1. However the fourth 
column of table 6 tell us this is not true. We have also studied this model 
with the twisted boundary conditions (0) and (§) with results as predicted 
in (|I0D, which clearly indicate an underlying c = 1 Gaussian field theory in 
the whole regime 1 > A > — 1. 



4.2 Model AB 2 

In this case, beyond the U(l) symmetry, we also have a Z(2) local gauge 
invariance corresponding to an independent interchange of the spin variables 
located at points like those shown in broken-lined rectangles in figure 2b. 
Since we have L/2 disjoint sectors labeled by the eigenvalues gi (±1) of the 
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gauge operators 

G l = a l .Bu + \ ; 1 = l,2,...,L/2, (26) 

where &i and are the spin-| operators located at the / th rectangle, a sector 
having gi = 1 will be spanned in a basis with three even combinations of the 
spin variables <?j and located at the retangle I. This means that in cr z -basis 

we should have the triplet combination | + +), ^=(| H — ) + | — h)) and | ). 

On the other hand if gi = — 1 we should have a singlet combination ^(| + 
— ) — | — h)). It is not difficult to verify that the interaction between spins 
in the sublattice A with a given neighbouring retangle I with gi = 1 (triplet) 
is exactly the same as in the (S, S') = (|, 1) Heisenberg interaction ( see ([!]) 
). In contrast, the interaction with a retangle with gi = —1 (singlet) is zero. 
This implies some interesting consequences. For a given U(l) sector S z = n 
the eigenenergies of the gauge sector with gi = 1 for all / = 1,2, . . . ,L/2 
of the Hamiltonian AB 2 with periodic ends will be exactly the same as the 
S z = n sector of the (S, S') = (^, 1) Heisenberg chain (|I|), also with periodic 
boundary condition. For general gauge choices we lose translation invariance 
since the operator (|5S|) does not have this symmetry. However, for the gauge 
choices g\ — g 2 — . . . — g^n — ±1) this invariance is recovered and we obtain 
the same dimensions as in the (S, S') = (|, 1) Heisenberg model studied in 
section 3. 

The eigenenergies in the sectors with g\ = —1 of the AB 2 model with 
periodic ends will correspond to the composition of energies of the (|, 1)- 
Heisenberg chains with free boundary conditions and different lattice sizes. 
This exact correspondence together with the relation ( |T2"D imply that the 
lowest energy in these sectors, in the bulk limit, will have a finite gap when 
compared with the ground-state energy, which happens in the sector g\ = 
gi — . . . — gi/2 — 1- This gap is proportional to the surface energy of 
the related (S,S') = (|, 1) Heisenberg chain. This produces the interesting 
feature that the correlation functions of operators which commute with ([2~6|) 
will have a power-law decay with exponents like those of the periodic (|, 1)- 
Heisenberg chain, while correlations of non-commuting operators may exhibit 
an exponential decay, with rate proportional to f^. An example of such 
operators is o%a^ af where k and / are indices inside a given retangle in Fig. 
2b and a = a x ± ia y . Apart from these pathological correlations most of 
them will be of the same nature as those of the (5*, 5") = (|, 1) Heisenberg 
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model and our results of section 3 indicate that they are described by a 
Gaussian like field theory in the regime 1 > A > — 1. 

5 Conclusions 

Anisotropic quantum chains with one kind of spin S exhibit a critical phase 
with continuously varying exponents governed by a c = 1 Gaussian like 
conformal field theory. This phase starts at the ferromagnetic point A = 1 
with an endpoint at A = A c (5), which for half-odd-integer S is expected 
to be 1 and A C (S) > —1 otherwise. This means that the antiferromagnetic 
point has quite different physics depending on the parity of 25. In this paper 
we analyse anisotropic Heisenberg chains with two kind of exchange-coupled 
centers. Due to a noncompensation effect, these models show ferrimagnetism 
instead of antiferromagnetism. As we change the anisotropy we move from 
the ferromagnetic (A = 1) to the ferrimagnetic (A = —1) point. We studied, 
by finite-size calculations and conformal invariance, four models of this kind; 
the (S, 5")-mixed Heisenberg chains with (S, S') = (|, 1) and (S, S') = (|, |), 
given in equation ([[J) and the Heisenberg models with XXZ interactions in 
the lattices of figures 2a and 2b ( models ABC and AB 2 ). We calculated the 
bulk and surface exponents of the first two models and the bulk exponents 
of the last two. 

All the models we studied show the universal feature of having a criti- 
cal phase for 1 > A > — 1 with long-distance physics governed by a c = 1 
Gaussian-like conformal field theory. The critical exponents, along this phase, 
are model-dependent continuous functions of the anisotropy. The sound ve- 
locity and compactification radius of the Gaussian theory go to zero at the 
isotropic ferromagnetic (A = 1) and ferrimagnetic (A = —1) points. This 
reflects the fact that at both points we should expect the appearance of 
quadratic dispersion relations. We strongly believe that this is the general 
scenario for arbitrary Heisenberg chains showing ferrimagnetism instead of 
antiferromagnetism. 
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Figure Captions 



Figure 1 - Schematic values of the lowest eigenenergy in a sector with mag- 
netization 5* 2 of the Hamiltonian ([]]). The ground-state energy is E and 
(S — S')^ is an integer. 

Figure 2 - Lattices where the (a) ABC and (b) AB 2 quantum chains are 
defined. At the circles we have spin-| SU(2) operators and along the lines 
the interactions are given by (^Sj) ( XXZ type ). The Hamiltonian AB 2 is in- 
variant under a local gauge transformation which independently interchanges 
the spin operators inside a rectangle in (b). 
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Table Captions 



Table 1 - Estimates, for some values of the anisotropy A, of the conformal 
anomaly of the (S, 5")-Heisenberg chain ([[]) for (S, S') = (~, 1) and (S, S') = 

a *) ~ 

V2' 2> 



Table 2 - Extrapolated values of the finite-size sequences ( [Eg) for the Hamil- 
tonian ([[J) with (S,S') = (|, 1) and some values of the anisotropy A. The 
extrapolations in the third and fourth column are obtained from the se- 

quences a;f' (L)/xi' (L) and (L)/x 2 ' 2 (L), respectively. The values of 

Xi' in the fifth column are obtained from standard finite-size scaling (see 
equation (|18|)). 

Table 3 - Extrapolated values of the finite-size sequences ( |T5| ) for the Hamil- 
tonian (|1|) with (S, S') = (|, |) and some values of the anisotropy A. The ex- 
trapolations in the third and fourth column are obtained from the sequences 

13 13 13 11 

x~2 ' 2 (L)/xi ' 2 (L) and xl" 1 (L) / xl" 2 (L) , respectively. 

Table 4 - Dimensions xf' s ($) obtained from the bulk limit of the finite- 
sequences (p0| ) obtained from the Hamiltonian (JJ) with twisted boundary 
conditions (|7],|[) specifyed by the angle $ = n and $ = The values 
in parentheses are the predicted ones obtained by using the values of xf' s 
estimated in tables 2 and 3 in (pip. 

Table 5 - Surface critical exponents xf s '(n) associated to the lowest eigener- 
gies in the sectors n = 1, 2 of the Hamiltonian ([!]) with (S, S') = (|, 1) and 
(S, S') = (|, |) for some values of A. These estimates are obtained from the 
sequences 

Table 6 - Anomalous dimensions x^ BC associated to the lowest eigenenergies 
in the sector n — 1, 2 of the model ABC defined in figure 2b, for some values 
of the anisotropy A. The third and fourth column are calculated similarly 
as for tables 2 and 3. 
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Table 1 



(S, S>) 


A = .8090 


A = 0.5 


A = 


A = -0.5 


A = -0.7071 


(1,1) 

\2' 2> 


0.93±0.05 
1.00±0.01 


1.01±0.01 
0.99±0.01 


1.03±0.04 
1.000±0.005 


1.01±0.01 
1.00±0.01 


0.9±0.1 
0.9±0.1 



Table 2 



A 


x( 


3.1 

x 2 


±1 ii 

xf ' jxl 


I,i II 

x \ / x i' 2 


1 - l P /2u 


0.9172 


0.0414±0.0001 


0.165±0.001 


3.9996±0.0005 


0.664±0.005 


0.038 


0.80901 


0.0639±0.0005 


0.255±0.001 


3.9999±0.0005 


0.653±0.008 


0.069 


0.5 


0.1041±0.002 


0.415±0.001 


3.999±0.005 


0.626±0.001 


0.104 


0.1736 


0.1302±0.0005 


0.521±0.004 


3.99±0.01 


0.584±0.002 


0.129 





0.139±0.001 


0.554±0.002 


3.999±0.004 


0.554±0.002 


0.138 


-0.5 


0.136±0.001 


0.547±0.005 


4.01±0.01 


0.40±0.01 


0.134 


-0.7071 


0.105±0.005 


0.44±0.02 


4.01±0.01 


0.31±0.01 


0.111 


-0.9010 


0.035±0.005 


0.15±0.02 


3.99±0.02 


0.11±0.01 


0.066 
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Table 3 





i 3 


1 3 


13 13 


13 11 


A 


x 1 


™2>2 


/■>-* 2 2/ rv> 2 2 
2 / 1 


x l 2 l x l 2 


0.9172 


0.0310±0.0005 


0.124±0.002 


3.9995±0.0005 


0.50±0.02 


0.80901 


0.047±0.001 


0.189±0.002 


4.001±0.001 


0.47±0.02 


0.5 


0.076±0.002 


0.301±0.002 


4.0000±0.0001 


0.44±0.01 


0.1736 


0.0914±0.0001 


0.3661±0.0002 


4.0001±0.0001 


0.412±0.001 





0.0959±0.0001 


0.384±0.001 


4.0000±0.0001 


0.383±0.001 


-0.5 


0.089±0.002 


0.357±0.003 


4.00±0.04 


0.26±0.01 


-0.7071 


0.071±0.002 


0.28±0.01 


3.98±0.01 


0.19±0.01 


-0.9010 


0.031±0.001 


O.lliO.Ol 


3.99±0.02 


0.065±0.008 



Table 4 



A = 0.8090 



A = 0.1736 



A = 



A = -0.5 



A = -0.7071 



xf'\<S> = 7T) 



Xi 2 ' ($ = 4tt/3) 



xH($ = 7r) 
= 4tt/3) 



0.93±0.02 
(0.978) 

1.65±0.05 
(1.74) 

1.26±0.05 
(1.33) 

2.22±0.05 
(2.36) 



0.49±0.01 
(0.480) 

0.86±0.01 
(0.853) 

0.686±0.002 
(0.684) 

1.222±0.004 
(1.216) 



0.46±0.01 
(0.4500 

0.81±0.01 
(0.799) 

0.657±0.005 
(0.652) 

1.165±0.005 
(1.159) 



0.48±0.01 
(0.460) 

0.84±0.02 
(0.817) 

0.70±0.01 
(0.702) 

1.25±0.01 
(1.248) 



0.58±0.01 
(0.595) 

1.03±0.02 
(1.06) 

0.90±0.03 
(0.880) 

1.53±0.04 
(1.565) 
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Table 5 



A 


x{ 


3.1 


1 3 


I 3 

™2>2 


0.9172 


0.082±0.002 


0.330±0.005 


0.061±0.001 


0.245±0.005 


0.80901 


0.128±0.001 


0.513±0.003 


0.098±0.003 


0.392±0.005 


0.1736 


0.258±0.002 


1.036±0.003 


0.183±0.001 


0.73±0.01 





0.272±0.005 


1.106±0.003 


0.191±0.001 


0.766±0.005 


-0.5 


0.266±0.002 


1.070±0.005 


0.181±0.002 


0.728±0.005 


-0.9010 


0.1086±0.0005 


0.440±0.002 


0.065±0.005 


0.260±0.002 



Table 6 



A 


^ABC 
x l 


^ABC 


„ABC I^ABC 
X 2 /X 1 


i i 

™ABC /~2'2 
x l 1 x l 


0.9172 


0.07±0.01 


0.27±0.02 


3.99±0.02 


0.86±0.03 


0.80901 


0.077±0.003 


0.31±0.01 


4.01±0.01 


0.84±0.02 


0.5 


0.129±0.002 


0.51±0.01 


4.00±0.03 


0.78±0.02 


0.1736 


0.158±0.003 


0.635±0.005 


4.01±0.02 


0.71±0.01 





0.165±0.002 


0.667±0.005 


4.01±0.02 


0.66±0.01 


-0.5 


0.152±0.002 


0.61±0.01 


3.99±0.02 


0.45±0.01 


-0.7071 


0.121±0.002 


0.486±0.002 


4.00±0.02 


0.32±0.01 


-0.9010 


0.05±0.01 


0.20±0.01 


3.99±0.03 


0.11±0.01 
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